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Measurements  of  the  spectral  and  intensity  dependences  of  the  optically-induced  reversal  of  current 
collapse  in  a  GaN  metal-semiconductor  field-effect  transistor  (MESFET)  have  been  compared  to 
calculated  results.  The  model  assumes  a  net  transfer  of  charge  from  the  conducting  channel  to 
trapping  states  in  the  high-resistivity  region  of  the  device.  The  reversal,  a  light-induced  increase  in 
the  trap-limited  drain  current,  results  from  the  photoionization  of  trapped  carriers  and  their  return  to 
the  channel  under  the  influence  of  the  built-in  electric  field  associated  with  the  trapped  charge 
distribution.  For  a  MESFET  in  which  two  distinct  trapping  centers  have  been  spectrally  resolved, 
the  experimentally  measured  dependence  upon  light  intensity  was  fitted  using  this  model.  The  two 
traps  were  found  to  have  very  different  photoionization  cross-sections  but  comparable 
concentrations  (4  X  1011  cm-2  and  6  X  1011  cm-2),  suggesting  that  both  traps  contribute  comparably 
to  the  observed  current  collapse.  [S0021-8979(00)00417-5] 


I.  INTRODUCTION 

Electronic  devices  based  on  GaN  and  AlGaN/GaN  het¬ 
erostructures  are  of  great  current  interest  because  of  the  po¬ 
tential  for  producing  structures  that  operate  at  high  tempera¬ 
ture,  high  power,  high  frequency,  and  in  adverse 
environments.  Significant  advances  have  been  made  in  the 
fabrication  of  these  devices:  Nitride -based  field-effect  tran¬ 
sistors  (FETs)  have  been  reported  to  exhibit  continuous  wave 
outputs  up  to  6.9  W/mm1  and  high  frequency  operation2  at 
fT=  67  GHz  and  /max=140GHz.  While  attainable  device 
characteristics  continue  to  improve,  the  reproducibility  of 
these  characteristics  has  become  problematic,3  due  to  device 
limitations  resulting  from  trapping  effects.3'4  In  addition  to 
compensating  intentionally  doped  material,  such  traps  can 
produce  persistent  photoconductivity  (PPC)  effects  and  cur¬ 
rent  collapse.  Current  collapse  refers  to  the  trapping  of 
charge  at  defects,  after  the  application  of  a  high  drain-source 
voltage,  which  results  in  a  significant  reduction  in  the  drain 
current.  This  effect  was  recently  investigated  in  a  GaN 
metal-semiconductor  field-effect  transistor5  (MESFET),  and 
was  observed  only  after  the  device  had  been  subjected  to  a 
large  drain  bias.  This  is  consistent  with  the  view  that  hot 
carriers  were  injected  into  the  high-resistivity  (HR)  GaN 
layer,  where  they  were  trapped  by  deep  centers.  Due  to  the 
low  thermal  emission  rate  from  these  defects,  the  trapped 
carriers  remain  long  after  the  applied  voltage  is  removed. 
This  transfer  of  charge  produces  a  depletion  region  in  the 
channel  of  the  device,  which  tends  to  pinch  off  the  device 
and  reduce  the  drain  current.  Current  collapse  has  been  a 
significant  problem  in  Si  metal-oxide-semiconductor  FETs,6 
in  AlGaAs/GaAs  modulation-doped  FETs,7  and  more  re¬ 
cently  in  nitride -based  AlGaN/GaN  heterostructure  FETs8 
(HFETs)  and  GaN  MESFETs.5’9 
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An  important  characteristic  of  current  collapse  is  that  it 
is  reversible  with  the  application  of  light:  the  trapped  carriers 
can  be  photoionized  and  are  thus  released  from  the  traps. 
Under  the  influence  of  the  built-in  electric  field  that  develops 
in  the  space-charge  region,  the  photoionized  carriers  drift 
rapidly  back  to  the  conducting  channel.  This  reduces  the 
width  of  the  depletion  region  in  the  channel  and  restores  the 
drain  current.  There  have  been  several  investigations  of  the 
influence  of  light  on  current  collapse.  Measurement  of  the 
spectral  dependence  of  the  conductance  in  AlGaAs/GaAs 
structures,  made  before  and  after  the  application  of  a  large 
drain-source  voltage,  led  to  the  conclusion10  that  the  traps 
responsible  for  current  collapse  were  localized  in  the  Al- 
GaAs  layer.  The  responsible  trap  was  later  identified  as  the 
DX  center  in  the  AlGaAs.7  A  later  probe  of  AlGaN/GaN 
devices  reported  increases  in  the  collapsed  drain  current  due 
to  illumination  with  near-band  gap  light  and  with  light  near 
650  nm,8  which  was  associated  with  an  unidentified  trap.  In 
studies  of  current  collapse  in  a  GaN  MESFET,  it  was 
observed5  that  the  efficiency  of  restoring  the  drain  current  by 
light  illumination  becomes  much  less  effective  with  increas¬ 
ing  wavelength.  While  current  interest  is  on  heterostructure 
devices  because  of  the  high  mobility  channel  formed  in  these 
structures,  it  is  also  important  to  understand  the  nature  of  the 
traps  in  the  underlying  GaN  in  order  to  be  able  to  unravel  the 
contributions  to  the  measured  optical  response  from  the  vari¬ 
ous  parts  of  the  device  structure.  Studying  current  collapse  in 
the  simpler  GaN  MESFET  provides  a  direct  way  of  carrying 
out  this  type  of  study. 

In  a  recent  report,9  we  measured  the  wavelength  depen¬ 
dence  of  the  optical  restoration  of  the  drain  current  in  a  GaN 
MESFET.  Under  the  initial  conditions  of  full  current  col¬ 
lapse  (i.e.,  before  illumination  and  with  all  the  traps  filled), 
the  fractional  increase  in  the  drain  current  per  incident  pho¬ 
ton  was  measured  as  a  function  of  wavelength.  This  optical 
response  function  is  related  to  the  photoionization  cross  sec¬ 
tion  of  the  trap.  The  resulting  spectrum  revealed  two  distinct 
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electron  traps.  Because  these  spectra  were  obtained  by  mea¬ 
suring  the  optical  reversal  of  current  collapse,  the  two  traps 
observed  must  be  those  responsible  for  this  effect.  Conse¬ 
quently,  this  type  of  spectroscopy  can  provide  signatures  for 
the  specific  traps  that  are  responsible  for  current  collapse  in 
FET  structures. 

In  this  work,  we  present  measurements  of  the  depen¬ 
dence  of  this  effect  upon  the  incident  light  intensity  and  de¬ 
velop  a  detailed  model  describing  the  optical  restoration  of 
the  collapsed  drain  current.  The  model  establishes  a  frame¬ 
work  for  understanding  the  previously  measured  spectral  de¬ 
pendence  and  for  analyzing  the  intensity-dependent  measure¬ 
ments.  This  analysis  also  provides  a  means  of  extracting 
important  trap  characteristics  from  the  experimental  data.  In 
addition,  the  model  determines  the  relationship  between  the 
spectral  dependence  of  the  measured  response  function  and 
that  of  the  trap  photoionization  cross  section. 

II.  EXPERIMENT 

A.  Sample  and  measurement  characteristics 

Details  of  the  MESFET  design  and  characterization  are 
described  in  Ref.  5.  The  FET  was  fabricated  with  a  drain- 
source  spacing  of  5  /rm,  a  gate  width  of  150  jam  and  a  gate 
length  of  1.5  jam.  The  200-nm- thick  n-GaN  active  channel 
layer  was  grown  on  top  of  a  3-/xm-thick,  undoped,  high- 
resistivity  GaN  buffer  layer  in  order  to  separate  the  active 
part  of  the  device  from  the  higher-defect-density  material 
near  the  sapphire  substrate.  Hall  measurements  at  300  K 
indicated  a  channel  carrier  concentration  of  2X  1 0 1 7 cm  3 
and  a  mobility  of  410  cm2/V s.  I—V  characteristics  were 
determined  using  a  HP4145b  semiconductor  parameter  ana¬ 
lyzer,  which  measures  the  drain  characteristics  with  a  single 
sweep  of  the  drain-source  voltage,  VD,  for  each  gate  bias. 
The  current  collapse  measurements  were  all  carried  out  at 
293  K  and  at  zero  gate-source  bias.  Monochromatic  light 
was  provided  by  a  75  W  Xe  arc  lamp  (PTI  model  A1010)  or 
a  tungsten-halogen  lamp,  and  a  Spex  1680B  0.22  m  double 
monochromator  with  a  1200  groove/mm  holographic  grating. 
The  spectrometer  bandpass  was  set  to  approximately  3.5  nm. 
The  light  was  collected  with  a  spherical  mirror  and  focused 
onto  the  device.  The  dependence  of  the  light-induced  in¬ 
crease  in  drain  current  upon  light  intensity  was  determined  at 
a  fixed  drain  voltage,  using  calibrated  neutral  density  filters 
to  attenuate  the  incident  light. 

B.  Measurement  of  drain  current  restoration 

Current  collapse  occurs  after  a  GaN  MESFET  is  sub¬ 
jected  to  a  high  drain-source  voltage,  and  hot  carriers  are 
injected  into  the  thick  HR  GaN  layer.  The  electrons  (for  an 
n -channel  device)  that  are  trapped  in  this  region  by  electri¬ 
cally  active  deep  centers  produce  a  charged  region  in  the  HR 
layer,  and  a  corresponding  depletion  layer  extends  into  the 
active  layer,  thus  narrowing  the  channel  and  causing  the  /-  V 
characteristic  of  the  device  to  collapse.  In  a  simplified 
“abrupt-junction”  view  of  this  process,  the  traps,  distributed 
uniformly  with  concentration  Nr,  fill  to  a  depth  xT  below  the 
n-type/HR  interface.  This  corresponds  to  the  case  in  which 
the  drain  current  is  fully  collapsed,  as  all  of  the  traps  within 


the  charged  region  are  filled.  As  the  device  is  irradiated  with 
light,  trapped  carriers  are  photoionized  and  drift  rapidly  back 
to  the  channel  under  the  influence  of  the  large  built-in  elec¬ 
tric  field  (s“105  V/cm,  see  Sec.  IV  C)  that  develops  across 
the  interface  due  to  the  charge  transfer.  This  reduces  the 
amount  of  trapped  charge,  shrinks  the  depletion  layer,  in¬ 
creases  the  drain  current,  and  reverses  the  current  collapse. 

As  the  fully  collapsed  (dark)  state  corresponds  to  the 
traps  being  completely  filled  within  the  charged  region  of  the 
HR-GaN,  a  natural  measure  of  the  reversal  of  this  phenom¬ 
enon  is  the  fractional  increase  in  the  drain  current  above  the 
dark  level.  This  is  a  direct  measure  of  the  fraction  of  trapped 
carriers  that  have  been  released.  To  be  meaningful,  this 
quantity  must  be  normalized  by  the  total  number  of  photons 
incident  on  the  sample.  Consequently,  at  a  fixed  value  of  the 
drain  voltage,  we  define  a  response  function9  that  is  the  frac¬ 
tional  increase  in  drain  current  above  the  dark  level  per  in¬ 
cident  photon  of  energy  hv 


S(hv)  = 


1 


Al(hv) 

^dark 


(1) 


Here  <I> hv  is  the  incident  photon  flux,  t  the  illumination  time, 
A I(hv)  the  light-induced  increase  in  the  drain  current  above 
the  dark  level;  Al(hv)  =  \I(hv)  —  /dark],  where  I{hv)  is  the 
drain  current  after  illumination  and  /dark  is  the  fully  collapsed 
drain  current.  The  product  •  t  is  the  total  exposure,  i.e., 
the  total  number  of  photons  per  unit  area  incident  on  the 
device  after  illumination  for  a  time  1.  In  Ref.  9,  we  noted  that 
under  certain  circumstances  the  response  function  S(hv)  is 
proportional  to  the  photoionization  cross  section  of  the  trap, 
cr (hv).  The  conditions  under  which  this  holds  is  determined 
in  Sec.  IV. 

In  order  to  ensure  reproducibility,  it  was  important  for 
the  measurements  to  begin  with  the  device  in  the  same  state 
for  each  measurement.  This  was  accomplished  by  irradiating 
the  device  with  a  blue  GaN  light  emitting  diode  (1  mW) 
before  applying  the  high  drain-source  bias.  This  procedure 
led  to  very  reproducible  results,  and  presumably  emptied  all 
of  the  traps,  thus  eliminating  the  channel  depletion  region 
and  returning  the  device  back  to  the  equilibrium  condition. 


C.  Photoionization  spectra 

In  order  to  provide  a  proper  context  for  the  intensity- 
dependent  measurements  and  for  the  model  that  will  be  de¬ 
veloped,  the  previously  reported9  spectral  dependence  of  the 
optical  response  function  S(liv)  is  shown  as  the  open  circles 
in  Fig.  1.  The  same  GaN  MESFET  is  used  in  the  present 
measurements.  Two  broad  absorptions  were  observed,  and 
labeled  Trap  1  and  Trap  2.  Standard  expressions  for  the  spec¬ 
tral  dependence  of  the  photoionization  cross  section  of  a 
deep  level11  were  unsuccessful  in  fitting  this  data.  From  the 
apparent  large  trap  depths  and  the  large  breadth  of  these 
absorptions,  it  was  concluded  that  the  traps  were  strongly 
lattice  coupled.  Consequently,  the  data  were  fitted  using  an 
expression  appropriate  for  absorption  from  strongly  lattice- 
coupled  defects,12  shown  as  the  solid  line  in  the  figure.  Two 
absorption  thresholds  were  determined,  at  1.80  and  2.85  eV 
below  the  conduction  band,  which  is  suggestive  of  deep 
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FIG.  1.  Spectral  dependence  of  the  293  K  response  function  S(hv),  which 
measures  the  fractional  increase  in  drain  current  above  the  dark  level  per 
incident  photon.  Contributions  to  the  absorption  are  evident  from  two  dis¬ 
tinct  traps.  The  solid  line  is  a  fit  using  the  spectral  dependence  of  the  photo¬ 
ionization  cross  section  appropriate  for  a  defect  coupled  strongly  to  the 
lattice  (Ref.  12).  The  dashed  lines  indicate  the  contributions  from  each  of 
the  two  traps. 

traps.  The  dashed  lines  represent  the  contributions  from  the 
individual  traps.  The  additional  rise  near  the  band  edge  was 
thought  to  be  due  to  the  optical  injection  of  carriers  into  the 
channel.  Photoionization  of  charged  shallow  acceptors  is 
also  a  possible  explanation.  Due  to  the  nature  of  the  mea¬ 
surement,  it  is  clear  that  the  two  traps  evident  in  Fig.  1  must 
be  those  responsible  for  current  collapse  in  the  device.  It  was 
also  noted9  that  these  same  traps  are  most  probably  associ¬ 
ated  with  PPC  effects  as  well,  as  each  of  the  two  spectra  was 
found  to  reproduce  the  photoconductivity  spectrum  of  a  trap 
associated  with  PPC.13'14 

In  order  to  study  each  trap  individually  through  the  de¬ 
pendence  upon  light  intensity,  a  fixed  wavelength  was  cho¬ 
sen  for  each  trap.  A  wavelength  of  600  nm  (2.07  eV)  was 
chosen  for  Trap  1  (see  Fig.  1).  At  this  wavelength,  the  photo- 
induced  increase  in  the  drain  current  is  still  reasonably  large, 
but  there  is  no  contribution  from  Trap  2.  The  spectrum  of 
Trap  2  overlaps  completely  with  contributions  from  Trap  1. 
However,  the  contribution  from  Trap  2  is  dominant  at  400 
nm  (3.10  eV).  Consequently,  this  wavelength  was  chosen  to 
characterize  Trap  2. 

III.  RESULTS 

Available  light  intensities  were  limited  to  ss  1 00  /xW 
spread  out  over  approximately  20  mm2.  The  illumination 
time  for  all  measurements  was  fixed  at  t  =  5  s  (as  in  Ref.  9), 
so  that  varying  the  light  intensity  is  equivalent  to  varying  the 
total  exposure,  (<FAuf).  In  what  follows  we  will  often  refer  to 
light  intensity,  the  independent  variable  in  the  experiment.  It 
should  be  kept  in  mind,  however,  that  the  change  in  the 
charge  density  due  to  the  optical  ejection  of  carriers  from  the 
traps  is  an  exposure  effect,  depending  on  both  the  light  in¬ 
tensity  and  the  illumination  time. 

As  the  dependence  of  A I{hv)  upon  light  intensity  was 
observed  at  a  fixed  wavelength,  it  was  equivalent  to  analyze 
either  A/(/u;)//dal.k  or  S(hv)  [see  Eq.  (1)],  as  long  as  the 


FIG.  2.  The  fractional  increase  in  the  drain  current  relative  to  the  dark  (fully 
collapsed)  level  is  plotted  (data  points)  against  the  incident  light  intensity  for 
600  nm  and  400  nm  excitation,  which  are  representative  of  Trap  1  and  Trap 
2,  respectively.  The  illumination  time  was  kept  constant  at  5  s.  The  solid 
lines  are  a  fit  to  the  data  using  Eq.  (22).  The  contribution  from  Trap  2  is 
represented  by  the  squares,  and  was  determined  by  subtracting  the  relatively 
small  contribution  of  Trap  1  at  400  nm  from  the  400  nm  Trap  2  data 
(crosses). 


duration  of  the  illumination  remained  fixed.  The  former  was 
chosen  as  a  more  intuitive  measure,  as  a  saturation  of 
A/// dark  is  expected  at  high  light  intensity  (exposure),  when 
almost  all  of  the  trapped  carriers  have  become  liberated.  In 
this  same  range  of  light  intensities,  S(hv )  decreases  as  (T>  1 
due  to  the  intensity  normalization  built  into  Eq.  (1).  The 
experimentally  measured  dependence  of  A///dark  on  the  inci¬ 
dent  intensity  is  shown  as  the  data  points  in  Fig.  2.  The  400 
nm  data  for  Trap  2  (crosses)  are  corrected  (squares)  by  sub¬ 
traction  of  the  small  contribution  from  Trap  1  at  400  nm,  as 
discussed  in  Sec.  IV.  A  linear  dependence  is  observed  at  low 
excitation  for  both  traps.  At  the  highest  excitations.  Trap  2 
appears  to  be  well  into  saturation,  while  the  response  of  Trap 
1  is  only  beginning  to  saturate.  This  is  consistent  with  the 
much  weaker  response  (by  a  factor  of  102)  of  Trap  1  com¬ 
pared  to  Trap  2  observed  in  the  spectral  study  shown  in  Fig. 
1.  As  we  will  see  in  Sec.  IV,  this  behavior  can  be  the  result 
of  either  a  smaller  photoionization  cross  section  or  a  lower 
trap  concentration  for  Trap  1. 

IV.  MODELING  OF  CURRENT  COLLAPSE 
A.  Single  trapping  species 

In  this  section,  we  will  model  current  collapse  in  the 
MESFET  in  order  to  develop  a  framework  for  understanding 
this  phenomenon  in  better  detail  and  to  enable  the  extraction 
of  physically  significant  parameters  from  the  experimental 
data.  We  reemphasize  that,  because  the  carriers  are  localized 
on  deep  traps,  the  spatial  distribution  of  charge  in  the  HR 
layer  is  effectively  “frozen  in”.  In  what  follows,  we  will 
emphasize  this  distinction  from  the  familiar  depletion  region 
by  using  the  term  ‘  ‘charged  region’  ’  instead.  This  distinction 
has  clear  physical  consequences:  in  a  standard  p-n  junction, 
as  the  amount  of  charge  in  the  depletion  layer  is  reduced  (for 
example,  by  applying  a  forward  bias),  the  width  of  the  deple¬ 
tion  layer  decreases  while  the  volume  charge  density  remains 
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FIG.  3.  A  sketch  of  the  GaN  MESFET  is  shown  with  the  parameters  rel¬ 
evant  to  the  model  calculation  in  Sec.  IV.  The  terms  h(y),  x„(y)  and  xT{y) 
represent  the  variations  along  the  channel  of  the  depths  of  the  gate  depletion 
width,  the  depletion  width  at  the  bottom  of  the  channel,  and  the  width  of  the 
charged  region,  respectively.  The  depth  of  the  active  region  is  d  and  the  gate 
length  is  L. 


the  same.  For  the  MESFET  considered  here,  this  remains  the 
case  only  for  the  positively  charged  depletion  layer  in  the 
active  region.  The  situation  for  the  charged  region  is  just 
reversed:  the  charge  within  the  region  is  trapped  at  deep 
defects.  If  some  of  these  trapped  carriers  are  liberated,  either 
optically  or  thermally,  and  return  to  the  channel,  the  remain¬ 
ing  charges  are  not  free  to  redistribute.  Consequently,  the 
layer  width  remains  the  same,  and  the  concentration  of 
trapped  charge  decreases. 

The  MESFET  will  be  modeled  in  a  manner  similar  to 
typical  textbook  treatments,15'16  except  that  these  assume 
symmetric  depletion  regions  associated  with  gate  contacts  at 
the  top  and  bottom  of  the  device.  We  will  assume  the  abrupt 
junction  approximation,  with  uniform  densities  of  dopants  in 
the  active  region  and  traps  in  the  HR  GaN.  The  latter  are 
assumed  to  be  neutral  when  empty  and  negatively  charged 
when  filled.  We  also  assume  the  gradual  channel  approxima¬ 
tion,  appropriate  for  a  1.5  /rm  gate  length,  and  with  the 
charge  distributed  uniformly  under  the  gate. 

The  device  geometry  is  depicted  in  Fig.  3.  Three  charged 
regions  exist:  There  are  depletion  regions  directly  under  the 
gate  and  at  the  bottom  of  the  channel  (due  to  the  trapped 
charge),  and  there  is  the  charged  region  in  the  HR  GaN  layer. 
After  the  application  of  a  large  drain-source  voltage,  the 
magnitude  of  the  initial  areal  charge  density  on  each  side  of 
the  n-type/HR  interface,  before  any  light  illumination,  is 
given  by  charge  conservation 

QT=eNdxn  =  eNTxT,  (2) 

where  e  is  the  magnitude  of  the  electronic  charge,  Nd  and  N  T 
are  the  (uniform)  volume  concentrations  of  donors  in  the 
channel  and  traps  in  the  HR  region,  respectively,  and  xn  and 
xT  are  the  initial  widths  of  the  depletion  and  charged  region 
layers,  respectively.  QT  represents  the  total  charge  trans¬ 
ferred  due  to  the  application  of  the  high  drain-source  voltage. 
It  is  assumed  that  all  of  the  traps  in  the  charged  region  are 
initially  filled.  After  illumination  with  light  for  a  time  t,  some 
of  the  trapped  charge  is  released  and  returns  to  the  channel, 


thus  decreasing  the  trapped  charge  density  and  the  width  of 
the  depletion  region.  The  dependence  of  the  areal  density  of 
trapped  charge  on  illumination  time  is  given,  following  Eq. 
2,  as 

Q(t)  =  eNdxn(t)  =  eN~(t)xT.  (3) 

The  depletion  layer  width  x„(t )  decreases  from  its  initial 
value  xn  as  the  total  amount  of  trapped  charge  is  reduced  by 
the  light  illumination.  The  width  of  the  charged  region  re¬ 
mains  fixed,  while  the  concentration  of  negatively  charged 
traps,  N~(t),  decreases  from  its  initial  value  NT  as  the  light 
releases  charge  from  some  of  the  traps.  The  concentration  of 
uncharged  traps  is  denoted  by  N°(t),  so  that 

N~(t)+N°(t)=NT.  (4) 

With  this  notation,  we  consider  the  modeling  of  current  col¬ 
lapse  in  the  MESFET. 

With  a  drain  voltage  applied,  the  voltage  drop  across  an 
incremental  section  of  the  channel  of  length  dy  and  resis¬ 
tance  dR  is 


dy 

dV=IDdR=ID-  — - — — -  (5) 

NdefxnW-[d-xn(y)-h(y)] 

where  //„  is  the  mobility  of  the  channel,  V(y  )  is  the  variation 
of  the  applied  potential  along  the  length  of  the  device,  W  the 
gate  width,  d  the  total  thickness  of  the  active  layer,  and  h  the 
width  of  the  depletion  region  under  the  gate.  The  expression 
in  square  brackets  is  just  the  width  of  the  conducting  region 
of  the  channel,  and  Nde/Jtn  is  the  channel  conductivity. 

Solving  Poisson’s  equation  in  one  dimension  and  with 
zero  gate  voltage  (the  experimental  condition)  leads  to  an 
expression  for  the  backside  depletion  width  in  the 
channel15'16 


x„(y)  = 


l+Nd/NTi 


■[Uu+V(y)l 


(6) 


where  t/bi  is  the  built-in  potential  across  the  n-type/HR  in¬ 
terface  and  e  is  the  static  dielectric  constant  of  the  material. 
The  width  of  the  depletion  layer  under  the  gate  is  given 

by15'16 


h{y)  = 


(7) 


with  Vbi  the  built-in  potential  associated  with  the  gate  con¬ 
tact.  Integrating  Eq.  (5)  over  the  length  of  the  gate,  L,  with 
V(0)  =  0  and  V(L)  =  VD ,  the  I-V  relationship  becomes 


7d(Fd)  — G 0- 1  VD 


2e 


3  V  eNdd 


(.vbi+vDy,z-v} 


3/2 _ i/3/2 

bi 


+ 


xll+Nd/NT 


=  (( Ubi+VD)m-U 5?) 


■,  (8) 


where  Gq  =  N de fxnWd/L  is  the  conductance  of  the  channel. 
At  low  drain  voltages,  VD<Vbi,  U b, ,  and  Eq.  (8)  becomes 
linear 
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Id~G  qVd 


1- 


1 

yll+Nd/NT 


(9) 

In  the  linear  regime,  the  depletion  widths  [Eqs.  (6)  and  (7)] 
become  independent  of  the  applied  drain  voltage,  and  depend 
only  on  the  built-in  potential.  The  measurements  described  in 
the  previous  section  were  carried  out  in  or  near  the  linear 
region.  Equation  9  expresses  the  I-  V  relationship  in  the  lin¬ 
ear  regime  and  with  no  light  irradiation:  This  corresponds  to 
the  drain  current  in  the  fully  collapsed  (dark)  state,  7dark .  The 
current  collapse  is  embodied  in  the  built-in  potential  t/bi, 
which  is  a  function  of  the  total  areal  charge  density  QT  [Eq. 
(2)]  injected  into  the  HR  region.  When  the  device  is  illumi¬ 
nated,  some  of  the  charge  is  released  from  the  traps  and 
returns  to  the  channel.  With  a  reduced  amount  of  trapped 
charge  in  the  HR  layer,  t/bi  also  decreases.  Thus,  the  built-in 
potential  at  the  interface,  and  therefore  the  drain  current,  de¬ 
pends  upon  the  amount  of  light  incident  upon  the  device: 
Uhl=  Uhl((l>J).  where  $  and  t  are  the  photon  flux  density 
(photons/cm2/s)  and  duration  of  the  incident  light  illumina¬ 
tion,  respectively.  The  increase  in  the  drain  current  (at  fixed 
V D)  due  to  the  light  is  obtained  from  Eq.  (9) 


A/(®,t)  =  /D(<I>,t)-Idaik 


Gov°  ^eNdd2(l+Nd/NTj 

•[Vt/bi(0,0)-Vt/bi($T)].  (10) 

Before  constructing  S(hv )  from  Eqs.  (1)  and  (10),  it  will 
be  convenient  to  express  the  built-in  potential  Uhl  in  terms  of 
the  areal  density  of  trapped  charge,  Q(t)  =  eNdx„(t)  [Eq. 
(3)],  where  the  depletion  width  in  the  channel  varies  with 
illumination  time  in  response  to  the  decreasing  trapped 
charge  density.  Combining  Eq.  (3)  with  Eq.  (6)  for  the  deple¬ 
tion  layer  width,  and  noting  that  in  the  linear  regime  x„  is 
independent  of  the  applied  voltage,  the  built-in  potential  with 
light  illumination  becomes 


Ubi(<S>,t)  = 


1+Nd/NT 

2eeNd 


(ID 


while  in  the  dark  we  have 


1  +Nd/NT 

U*{0,0)  =  2eeNd  [8(0,0)]  ,  (12) 

where  Qj=Q{ 0,0)  corresponds  to  the  total  areal  charge  den¬ 
sity  initially  transferred  into  the  HR  region.  Using  Eqs.  (1) 
and  (9)-(12),  the  response  function  for  photon  flux  density 
at  photon  energy  hv  for  illumination  time  t  then  be¬ 
comes 


S(hv,t)  = 


1 

®hv-t 


Qt  Qi'&hv’t) 
eNdd—  \]2e  eN  dVhi—  QT 


(13) 


The  denominator  of  the  expression  in  brackets  reflects  the 
areal  density  of  channel  carriers,  after  the  depletion  regions 
under  the  gate  and  at  the  HR  interface  are  taken  into  account. 

A  determination  of  the  response  function  S(hv )  requires 
evaluating  the  dependence  of  the  trapped  charge  density 
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FIG.  4.  Schematic  level  diagram  indicating  the  dominant  processes  that 
determine  the  dynamics  of  the  photoexcited  carrier  population  in  the  HR 
GaN  region:  photoionization,  carrier  capture,  and  carrier  drift  out  of  the 
charged  region.  The  illumination-time  dependent  carrier  concentration,  neu¬ 
tral  trap  and  charged  trap  concentrations  are  given  by  n(t),  and 

respectively.  Terms  describing  the  transition  rates  for  each  process 
are  also  shown. 


upon  the  illumination  time  and  intensity.  At  t  =  0,  just  before 
illumination  begins,  the  traps  are  filled  within  the  charged 
region,  N~(0)=NT.  The  areal  density  of  trapped  negative 
charge  after  light  illumination  for  duration  t  is  [Eq.  (3)] 
Q(t)  =  eN~(t)xT,  where,  as  we  noted  earlier,  xT  remains 
fixed  because  the  remaining  trapped  charge  cannot  redistrib¬ 
ute.  The  total  charge  in  the  HR  layer  should  also  include 
contributions  from  the  photoexcited  carriers  in  the  HR  region 
that  have  not  yet  crossed  back  into  the  channel.  This  contri¬ 
bution  will  be  neglected  because,  as  we  will  show  in  Sec.  V, 
this  concentration  is  extremely  low,  as  the  large  built-in  field 
drives  the  photoexcited  carriers  rapidly  back  into  the  channel 
layer. 

In  equilibrium,  (no  trapped  charge)  the  traps  can  be  ei¬ 
ther  empty  or  partially  compensated,  depending  upon  the  po¬ 
sition  of  the  Fermi  level  in  the  HR  material.  If  the  traps  are 
partially  compensated,  the  total  trap  concentration  NT  in  Eq. 
(2)  must  correspond  only  to  the  concentration  of  uncompen¬ 
sated  traps  in  equilibrium.  The  remaining  traps  that  were 
already  filled  in  equilibrium  will  not  contribute  to  the  current 
collapse  process  or  to  the  optical  reversal  of  current  collapse, 
since  these  processes  only  depend  upon  the  amount  of  charge 
that  has  been  transferred  from  the  channel  to  the  HR  layer. 
For  example,  if  all  NT  trapped  carriers  have  been  photoion- 
ized  and  removed  from  the  charged  region,  any  remaining 
filled  traps,  corresponding  to  those  that  were  already  filled  in 
equilibrium,  can  still  be  photoionized.  However,  these  pho¬ 
toexcited  carriers  experience  no  built-in  field,  as  the  region  is 
no  longer  charged.  Consequently,  they  remain  in  the  HR 
region  and  are  eventually  re -trapped.  With  these  consider¬ 
ations,  N~{t)  must  correspond  only  to  the  traps  that  were 
filled  as  a  consequence  of  the  initial  application  of  a  high 
drain-source  voltage. 

The  carrier  dynamics  associated  with  the  photoionization 
of  traps  in  the  charged  region  is  represented  in  Fig.  4.  Three 
processes  dominate  the  time  dependence  of  the  trapped 
charge  density:  photoionization,  carrier  capture,  and  carrier 
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drift  back  to  the  conducting  channel.  The  former  process 
generates  photoexcited  carriers  into  the  conduction  band  of 
the  HR  material,  while  the  latter  two  processes  compete  for 
those  carriers.  In  Sec.  V,  we  will  show  that  for  the  conditions 
under  which  the  experiments  were  carried  out,  the  carrier 
capture  rate  is  negligible  in  comparison  to  the  rate  of  carrier 
drift.  Neglecting  carrier  capture,  N~(t)  takes  a  particularly 
simple  form 

N-{t)  =  NTe-‘T(hu]'t,\  (14) 

where  cr(hv)  is  the  photoionization  cross  section  of  the  trap. 
Consequently,  we  can  rewrite  Eq.  (3)  as 

Q{^,t)  =  eNTxTe~,r(hv',mt=QTe-,T(hv)‘s,t.  (15) 

Combining  Eqs.  (13)  and  (15),  the  dependence  of  the  re¬ 
sponse  function  on  illumination  time  and  intensity  becomes 

S(hv,t)=^--K-[  i-g-^-®"]  (16) 

hv  t 

with 


K= 


Qt 

eNdd~  'j2eeNdVbi~QT 


(17) 


where  the  parameter  K  is  dependent  only  upon  the  initial 
trapped  charge  density  and  the  characteristics  of  the  channel. 
Under  conditions  of  weak  carrier  excitation  (small  photoion¬ 
ization  cross  section  or  low  light  exposure),  such  that 
cr{hv)  ■  •  f<  1,  Eq.  16  becomes 

S{hv)=K-  cr(hv).  (18) 


Thus,  the  conditions  under  which  the  response  function 
S(hv)  is  proportional  to  the  photoionization  cross  section  of 
the  trap  are;  (1)  that  the  measurement  is  made  in  or  near  the 
linear  regime,  VD<Vhi,  Ubl,  and  (2)  that  the  measurement  is 
made  under  weak  excitation,  cr{hv)  •  <3?  •  1.  The  photoion¬ 

ization  spectra  shown  in  Fig.  1  were  obtained  under  these 
conditions. 

It  is  now  evident  that  the  response  function  S(hv )  is  not 
sufficient  to  determine  whether  one  of  the  two  traps  in  Fig.  1 
is  dominant  in  causing  the  observed  current  collapse.  While 
the  spectra  in  Fig.  1  clearly  resolve  the  responsible  traps,  the 
response  function  only  measures  how  effectively  light  can 
reverse  this  effect.  The  parameter  that  determines  the  contri¬ 
bution  of  a  particular  trap  toward  producing  current  collapse 
is  the  concentration  of  carriers  that  can  be  trapped  on  that 
defect  due  to  the  application  of  the  high  drain-source  volt¬ 
age.  Equations  (2),  (17)  and  (18)  indicate  that  S(hv )  depends 
on  both  the  trap  concentration  and  on  the  photoionization 
cross  section.  These  parameters  are  deduced  in  Sec.  IV  C  by 
fitting  the  experimental  data  to  the  model  developed  here. 


B.  Multiple  trap  species 

For  wavelengths  longer  than  about  480  nm  ( hv 
<2.6  eV),  the  contribution  of  Trap  1  to  S(hv )  in  Fig.  1  is 
independent  of  that  from  Trap  2.  For  shorter  wavelengths, 
both  traps  contribute.  The  model  described  above  can  be 
extended  rather  simply  to  the  case  of  overlapping  contribu¬ 
tions,  such  as  that  for  liv> 2.6eV  in  Fig.  1.  Assuming  that 


the  charge  is  distributed  over  two  distinct  traps,  it  is  straight¬ 
forward  to  show  that  the  response  function  becomes 

S(hv,t)=^--{KV[  l-g-^fcO-*-*] 

hv  t 

+  Kr[  (19) 

where 

Kj  = -  QjT  - .  (20) 

eNdd--j2eeNdVh-Qr 

Here,  the  (Tj(hv)  and  the  QjT  are  the  photoionization  cross 
section  and  the  areal  charge  density  associated  with  each 
trap,  y  =  l,2,  and  Qt=  Q\t+  Qit  is  the  total  areal  charge 
density  [Eq.  (2)].  Extension  to  more  than  two  traps  is  obvi¬ 
ous.  For  the  case  of  low  excitation,  crj(lw )  •  <b  -  f  <  1 ,  and  the 
response  function  reduces  to 

S(hv)  =  Kl  •  crl(hv)  +  K2-  cr2(hv).  (21) 

Thus,  the  individual  contributions  in  Fig.  1  are  additive  un¬ 
der  low  excitation  conditions. 


C.  Fitting  of  the  intensity  dependence 

The  dependence  of  the  optical  reversal  of  current  col¬ 
lapse  on  the  intensity  of  the  incident  light  for  fixed  illumina¬ 
tion  time  (Fig.  2),  was  studied  by  monitoring  A/(/u;)//dark 
rather  than  S(hv),  as  noted  in  Sec.  III.  From  Eqs.  (1)  and 
(16),  this  can  be  written 


A  I(hv) 


1  dark 


=*•[1 


e-o-(/iu)$r 


]. 


(22) 


where  K  is  a  function  of  the  concentration  of  trapped  charge 
and  known  channel  characteristics  [Eq.  (17)].  At  400  nm,  the 
response  seen  in  Fig.  2  for  Trap  2  (crosses)  contains  a  small 
contribution  from  Trap  1  (see  Fig.  1)  that  becomes  non- 
negligible  in  the  saturation  regime.  Subtracting  this  contribu¬ 
tion,  the  squares  in  Fig.  2  indicate  the  response  associated 
with  Trap  2  alone.  For  each  trap,  two  parameters  are  required 
to  fit  the  data;  K  and  cr(hvexp),  where  h  v  exp  is  the  photon 
energy  used  in  the  measurement.  Equation  (22)  is  plotted  as 
the  solid  lines  in  Fig.  2,  using  the  results  of  this  fitting  pro¬ 
cedure  for  each  of  the  two  traps.  The  parameter  K  determines 
the  saturation  level,  corresponding  to  having  emptied  all  of 
the  traps,  while  the  magnitude  of  the  data  in  the  linear  region 
is  most  sensitive  to  the  photoionization  cross  section.  The 
resulting  parameters  were  found  to  be:  cr2(400  nm) 
s2.8Xl(T15  cm2,  K2=0A0 ;  o-,(600  nm)s6X10-17  cm2, 
AT !  =  0.24.  Note  that  the  relative  photoionization  cross  sec¬ 
tions  for  the  two  traps  as  determined  from  this  procedure  are 
consistent  with  the  spectral  measurements  of  Fig.  1. 

From  the  experimentally  determined  values  of  K  for 
each  trap,  the  areal  trap  concentrations  A flT  and  AT2t  (MjT 
=  QjTle=NjT-xT )  can  be  determined  using  Eq.  20.  With 
Nd=  2X1017  cm-3,  d  =  200  nm,  Vbi=l  V  for  the  Pt  gate 
contact,  and  taking  es9.5  for  GaN, 17-19  we  find  that  J\flT  and 
J\f2r  are  comparable:  J\flT=4X  1011  cm-2,  A/^T-sbXlO11 
cm-2,  and  consequently  the  total  areal  trap  density  is  J\fT 
six  10l2cm  2.  It  is  now  clear  that  the  large  difference  in 


J.  Appl.  Phys.,  Vol.  88,  No.  5,  1  September  2000 


Klein  et  at. 


2849 


the  observed  magnitudes  of  S(hv )  for  the  two  traps  (Fig.  1) 
results  primarily  from  the  difference  in  their  photoionization 
cross  sections.  In  contrast,  the  concentrations  of  the  two  traps 
are  comparable.  Consequently,  both  traps  contribute  compa¬ 
rably  to  the  current  collapse  in  this  device. 

The  maximum  field  across  the  n-type/HR  interface  can 
now  be  estimated:  for  the  simple  abrupt  junction  described 
earlier,  Smax=QTle=e-J\fTl  e.  We  find  the  field  at  the  inter¬ 
face,  £max=2X  105  V/cm  for  this  device  in  the  fully  collapsed 
state.  Thus,  the  transfer  of  charge  into  the  HR  region  pro¬ 
duces  a  rather  large  initial  electric  field  that  drives  the  carri¬ 
ers  rapidly  back  into  the  channel  when  they  are  released  from 
the  traps. 

V.  CARRIER  DYNAMICS 

Two  simplifying  assumptions  were  made  in  connection 
with  Eqs.  (3)  and  (14)  in  determining  the  dependence  upon 
illumination  time  of  the  concentration  of  filled  traps,  N~{t), 
and  of  the  total  areal  charge  density,  Q{t).  These  assump¬ 
tions  were:  (1)  that  the  rate  of  recapture  of  photoexcited  car¬ 
riers  in  the  HR  region  was  negligible  in  comparison  to  the 
rate  of  drift  of  carriers  out  of  the  region,  and  (2)  that  the 
charge  density  due  to  photoexcited  free  carriers  in  the  HR 
region  was  negligible  compared  with  the  charge  density  due 
to  trapped  carriers.  It  will  be  shown  that  these  assumptions 
are  valid  except  in  the  limit  of  very  long  illumination  times, 
when  almost  all  of  the  traps  have  been  emptied,  and  the 
built-in  field  has  become  weak.  If  the  carrier  generation  rate, 
<t{  h  v )  (F ,  becomes  small  enough  (low  optical  fluence  or 
small  photoionization  cross  section),  thermal  emission  would 
also  have  to  be  taken  into  account. 

To  examine  this  situation  more  closely,  the  dynamics  of 
carrier  generation  and  recombination  must  be  considered. 
The  photoexcited  carrier  dynamics  in  the  HR  GaN  was  sum¬ 
marized  in  Fig.  4.  Trapped  carriers,  with  concentration 
N~(t),  are  photoionized  by  the  incident  light,  and  may  be 
re-trapped  at  neutral  (uncompensated)  centers,  with  concen¬ 
tration  N°{t)  =  NT—N~{t),  or  they  may  be  swept  out  of  the 
charged  region  by  the  built-in  field  before  they  can  be 
trapped.  The  current  density  J„  associated  with  the  photoex¬ 
cited  carriers  can  be  related  to  the  carrier  concentration  n(t ) 
in  the  HR  region  through  the  continuity  equation 

dn  1 

—  =  Gn-Rn+-\Jn,  (23) 

dt  e 

where  Gn  is  the  carrier  generation  rate  due  to  photoioniza¬ 
tion  and  Rn  is  the  recombination  rate  due  to  re-trapping.  As 
indicated  in  Fig.  4,  we  may  write 

Gn  =  <r(hv)<bN~(t),  (24) 

R„  =  <rcVtoN0(t)n(t),  (25) 

where  crc  is  the  room  temperature  electron  capture  cross 

section  of  the  trap,  uth  is  the  average  electron  thermal  veloc- 

ity,  nth=(3 kBT/m*)m,  kB  is  Boltzmann’s  constant,  and  m* 
is  the  electron  effective  mass  of  GaN,20  m*  =  0.2m0.  Be¬ 
cause  of  the  high  electric  field  in  the  charged  region,  carrier 
drift  effects  will  dominate  carrier  diffusion.  Under  these  con¬ 
ditions,  Eq.  (23)  simplifies  to 


dn  d£ 

-  =  Gn-R„  +  np.e—,  (26) 


where  £  is  the  built-in  electric  field  and  pe  the  electron  mo¬ 
bility  in  the  HR  GaN.  Note  that  the  very  low  mobility  ob¬ 
served  in  HR  GaN,  generally  associated  with  a  deep-defect 
hopping  conductivity  mechanism,21  is  appropriate  only  for 
thermally  generated  carriers.  The  mobility  of  photoexcited 
carriers  would  be  expected  to  be  closer  to  that  measured  in 
the  n-type  layer. 

Under  the  abrupt  junction  approximation  with  uniform 
charge  densities,  £{x)  is  linear  in  x.  Taking  the  total  volume 
charge  density  in  the  HR  region  as  —  e\N~{t)  +  n{t)], 
which  allows  contributions  from  both  charged  traps  and  pho¬ 
toexcited  carriers,  we  find  that  d£ldx=  —  e\N~(t) 
+  n(t)\le.  Combining  this  with  Eqs.  (24) — (26),  the  time- 
dependent  carrier  concentration  is  determined  by 

dn 

—  =  a(hv)&N  ( t)  —  {aN°(t)  +  b[N  ( t)  +  n(t)]}n(t ), 

(27) 

where  a  =  crcu  *  and  b={  peel  e)  are  representative  of  the 
capture  and  drift  processes,  respectively.  The  competition 
between  capture  and  drift  is  evident  in  the  term  enclosed  in 
curly  brackets.  The  corresponding  concentration  of  nega¬ 
tively  charged  traps  is  determined  by 


dN~ 

dt 


—  cr{hv)^N  ( t)  +  aN°(t)n(t ). 


(28) 


If  we  first  consider  the  case  where  recapture  of  the  pho¬ 
toexcited  carriers  is  neglected  (a  =  0),  Eq.  (28)  leads  directly 
to  the  exponential  decay  of  the  trapped  carrier  concentration 
indicated  in  Eq.  (14).  The  solution  of  Eq.  (27)  for  n(t)  in  this 
situation  results  in  a  carrier  concentration  that  rises  rapidly 
with  a  time  constant  (bNT)~x  to  a  constant  level,  nss 
=  [o-{hv)‘t>/b],  determined  by  a  quasisteady-state  balance 
between  carrier  generation  and  carrier  drift.22  Assuming  a 
typical  trap  concentration  of  1 0 1 6  cm  3  for  Trap  2 
[o-2(400 nm)«*2.8X  1CU15 cm2],  the  carrier  concentration  is 
found  to  rise  with  a  characteristic  time  of  about  1  psec  to  a 
maximum  quasisteady-state  value  (for  the  highest  optical  flu¬ 
ence  employed  in  this  work)  of  approximately 
nss®s6Xl04cm-3.  This  is  shown  as  the  inset  in  Fig.  5  for  an 
incident  photon  flux  of  <F  =  5X  1013cm_2s_1  (see,  e.g..  Fig. 
2).  It  thus  becomes  clear  from  the  magnitude  of  nss  that  the 
charge  density  associated  with  the  free  carrier  concentration 
remains  much  smaller  than  that  associated  with  the  filled 
traps,  until  the  latter  are  almost  all  empty.  Consequently,  the 
neglect  of  the  free  carriers  in  determining  the  charge  density 
of  the  HR  region  is  seen  to  be  a  very  reasonable  approxima¬ 
tion. 

For  carrier  capture  to  be  negligible,  the  capture  term  in 
Eq.  (27)  must  be  much  smaller  than  the  drift  term.  Taking 
into  account  that  n{t)<NT,  this  condition  becomes 


a\ 

N°(t) 

W 

[at(o] 

The  condition  under  which  this  inequality  remains  valid  is 
determined  by  the  magnitude  of  {alb).  Evaluating  this  ex- 
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FIG.  5.  The  dependence  upon  illumination  time  of  the  photoexcited  carrier  concentration,  n(t ),  and  the  charged  trap  concentration,  N  ( f ) ,  calculated  by 
solving  Eqs.  (27)  and  (28)  simultaneously  by  numerical  methods.  The  top  row  (a,b,c)  corresponds  to  a  trap  with  a  negligible  capture  cross  section,  while  the 
bottom  row  (d,e,f)  corresponds  to  a  room  temperature  capture  cross  section  crc=  1  X  10“ 18  cm2.  Three  values  of  the  incident  photon  flux  <I>  are  indicated, 
corresponding  with  increasing  excitation  to  linear,  intermediate,  and  saturation  regimes,  as  observed  for  Trap  2  in  Fig.  2.  The  optical  exposure  in  these 
calculations  is  varied  at  fixed  light  intensity  by  varying  the  illumination  time.  The  inset  shows  the  development  of  the  quasisteady-state  carrier  concentration 
on  the  picosecond  time  scale,  calculated  for  intermediate  excitation  conditions,  <I>  =  5  X  1013  cm~2  s_1. 


pression  requires  knowledge  of  the  room  temperature  capture 
cross  section  of  the  traps  being  investigated  in  this  work. 

While  such  data  are  not  available  for  the  traps  consid¬ 
ered  here,  capture  cross  sections  in  HR  and  n-type  GaN  have 
been  reported23  25  in  a  broad  range  from  roughly  10” 23  cm2 
to  10” 15  cm2.  Polyakov  et  al.25  have  reported  capture 
cross  sections  for  GaN  near  room  temperature  for  the  E2 
center  [trc(278 K)  =  4 X 10”20  cm2]  and  the  E4  center 
[crc(385  K)  =  3.3X  10”23  cm2].  In  Ref.  9  we  presented  evi¬ 
dence  to  suggest  that  the  two  traps  observed  in  this  work  also 
participate  in  PPC.  As  PPC  due  to  traps  generally  involves 
centers  with  small  capture  cross  sections,  it  follows  that  the 
traps  in  Fig.  1  would  also  be  expected  to  have  small  capture 
cross  sections.  Hirsch  et  al. 13  confirm  this,  concluding  from 
their  PPC  measurements  that  the  trap  corresponding  to  Trap 
2  in  the  current  work  would  need  to  have  a  very  small  cap¬ 
ture  cross  section.  With  these  considerations  in  mind,  we 
choose  a  capture  cross  section  of  1  X  10” 18  cm2  as  an  appro¬ 
priate  compromise  to  estimate  {alb)  for  Trap  2.  This  value  is 
small  compared  to  the  largest  capture  cross  sections,  yet  not 
so  small  as  to  preclude  any  significant  re-trapping. 

Using  this  value  for  crc ,  a  mobility  comparable  to  that  in 
the  n-GaN  channel,  and  with  the  300  °K  thermal  velocity  for 
GaN,  uths9xl06cm/s,  we  find  (a/b)^  1X10”7.  This  sug¬ 
gests  that,  even  with  uncertainty  in  some  of  the  parameters, 
the  inequality  given  by  Eq.  (29)  remains  valid  until  virtually 
all  of  the  charged  traps  have  been  neutralized.  Consequently, 
the  neglect  of  carrier  trapping  is  of  concern  only  at  very  long 
illumination  times  and  does  not  have  a  significant  effect  on 
the  analysis  of  the  data  in  Fig.  2. 


These  considerations  may  be  investigated  in  more  detail 
by  solving  the  coupled  equations  (26)  and  (27)  simulta¬ 
neously,  and  without  further  approximations,  using  numeri¬ 
cal  methods.  The  result  of  this  computation  is  shown  in  Fig. 
5.  The  concentrations  of  N~{t)  (dashed  line)  and  n(t)  (solid 
line)  are  plotted  versus  illumination  time  for  three  values 
of  the  incident  light  flux;  5  X  1014,  5  X  1013,  and  1  X  1013 
cm”2s_1,  which  correspond  to  saturation,  intermediate  and 
linear  regimes,  respectively,  for  Trap  2  in  the  intensity  de¬ 
pendent  measurements  shown  in  Fig.  2.  The  top  row  of  pan¬ 
els  corresponds  to  the  case  with  negligible  re-trapping,  while 
the  bottom  row  corresponds  to  a  room  temperature  capture 
cross  section  of  IX  10” 18  cm2.  For  the  time  scale  in  this 
figure,  the  carrier  concentration  is  initially  constant  at  its 
quasisteady-state  value,  \a{hv)‘t>/b],  as  shown  in  the  inset. 
Without  carrier  capture,  the  trapped  carrier  concentration  de¬ 
cays  exponentially  with  irradiation  time,  while  the  carrier 
concentration  remains  constant  until  the  filled  trap  concen¬ 
tration  has  decayed  to  the  point  where  N~(t)^n(t).  At  this 
stage,  N~(t)  no  longer  dominates  the  contribution  to  the 
charge  density  that  determines  the  built-in  electric  field.  Be¬ 
yond  this  point  the  quasisteady-state  carrier  concentration 
can  no  longer  be  supported,  and  n(t )  begins  to  decay.  As 
expected,  more  intense  light  leads  to  a  faster  decay  of  the 
filled  trap  concentration,  and  thus  to  an  earlier  decay  of  the 
carrier  concentration. 

The  bottom  row  of  panels  in  Fig.  5  corresponds  to  a 
nonzero  capture  cross  section.  There  is  little  effect  of  capture 
on  the  carrier  and  filled  trap  concentrations  at  short  illumi¬ 
nation  times.  Capture  effects  become  significant  at  later 
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times,  when  the  inequality  in  Eq.  (29)  is  no  longer  satisfied. 
It  can  be  shown  that  this  same  condition  also  corresponds  to 
N~(t)  becoming  nonexponential,  as  the  recapture  of  photo- 
ionized  carriers  slows  the  decay  of  the  trapped  carrier  con¬ 
centration.  With  stronger  excitation,  capture  becomes  signifi¬ 
cant  at  earlier  illumination  times.  At  later  times,  when  there 
is  significant  re-trapping  of  the  photoionized  carriers,  the 
quasisteady-state  balance  of  photoionization  and  drift  is  dis¬ 
rupted,  and  the  carrier  concentration  begins  to  decay.  This 
occurs  when  the  fraction  of  filled  traps,  N  (t)/Nr,  has  de¬ 
creased  sufficiently  so  as  to  be  comparable  with  the  fraction 
(alb),  and  almost  all  of  the  traps  have  been  emptied.  At 
even  later  times,  when  capture  dominates  over  drift,  n(t) 
follows  a  r1  behavior.  It  should  be  noted  that  these  consid¬ 
erations  apply  to  time  scales  that  are  substantially  longer 
than  the  5  s  illumination  time  employed  in  the  measurements 
of  Figs.  1  and  2.  The  experimental  conditions  correspond  to 
a  quasisteady-state  carrier  concentration  and  an  exponen¬ 
tially  decaying  filled  trap  concentration. 

VI.  DISCUSSION 

The  experimental  technique  discussed  in  this  work  has 
been  shown  to  provide  spectroscopic  signatures  of  traps  that 
degrade  FET  output  characteristics  by  limiting  the  drain  cur¬ 
rent.  By  modeling  the  current  collapse  process,  important 
trap  characteristics  such  as  the  areal  trap  density  and  the 
photoionization  cross  section  can  be  obtained.  Of  greater  sig¬ 
nificance,  the  model  provides  a  framework  for  understanding 
the  process  and  for  applying  the  technique.  The  main  limita¬ 
tion  of  the  measurement,  of  course,  is  that  it  is  sensitive  only 
to  current  collapse  and  to  traps  that  participate  in  that  pro¬ 
cess.  In  some  sense,  this  is  also  an  attribute,  as  it  is  clear  that 
any  traps  observed  in  this  way  can  be  unambiguously  iden¬ 
tified  as  contributors  to  this  effect.  This  limitation  should  not 
be  considered  absolute,  in  the  sense  that  the  general  approach 
may  be  extended  to  any  trap-related  phenomenon  that  re¬ 
sponds  to  light  illumination.  Photocapacitance  spectroscopy 
is  an  obvious  example. 

When  extracting  physical  parameters  using  the  model 
developed  in  Sec.  IV,  it  is  most  suitable  to  use  a  gateless 
structure,  as  the  gate  tends  to  block  some  of  the  light,  thus 
leading  to  an  overestimate  of  the  incident  light  intensity. 
When  such  structures  are  not  available,  as  in  the  present 
case,  it  can  be  shown  that  this  effect  will  lead  to  an  under¬ 
estimate  of  the  extracted  photoionization  cross  sections  by 
the  same  factor,  but  will  have  no  significant  effect  on  the 
extracted  areal  trap  concentrations  or  on  the  spectral  depen¬ 
dence  of  the  photoionization  cross  section  (Fig.  1). 

While  the  use  of  a  gateless  structure  would  be  ideal, 
there  are  several  approaches  that  could  be  employed  with 
gated  structures.  One  solution  is  to  use  a  semitransparent 
gate.  While  FETs  are  not  generally  fabricated  in  this  way,  it 
would  be  possible  to  produce  test  structures  in  this  manner. 
As  the  substrates  are  transparent  to  light,  another  alternative 
would  be  to  allow  the  light  to  enter  from  the  bottom  of  the 
device.  However,  this  limits  the  available  wavelengths  to 
those  below  the  band  gap  of  GaN,  which  could  become  prob¬ 
lematic  in  studies  of  AlGaN/GaN  structures.  Bottom-access 


measurements  have  been  used  to  probe  the  actual  location  of 
traps  in  a  device  structure:  Dang  et  al.26  measured  the  spec¬ 
tral  dependence  of  photoconductivity  in  a  PPC  study  of  an 
AlGaN/GaN  structure,  comparing  results  from  top  and  bot¬ 
tom  access  in  an  effort  to  determine  the  location  (AlGaN  or 
GaN)  of  traps  associated  with  features  in  the  photoconduc¬ 
tivity  spectrum.  It  was  concluded  that  defects  giving  rise  to 
PPC  were  present  in  both  the  AlGaN  and  the  GaN  layer. 
While  some  features  could  be  associated  with  a  particular 
layer,  the  location  of  other  spectral  features  that  were  below 
the  gap  of  both  materials  could  not  be  distinguished.  The 
broad  absorption  observed  at  longer  wavelengths  was  inter¬ 
preted  in  terms  of  a  broad  distribution  of  defect  states  within 
the  band  gap.  Such  a  broad  distribution  is  not  required,  how¬ 
ever,  in  order  to  account  for  this  kind  of  absorption  spec¬ 
trum.  As  we  have  noted9  in  reference  to  the  two  broad  ab¬ 
sorptions  in  Fig.  1,  such  broad  bands  can  result  quite 
naturally  from  a  single  deep  defect  that  is  strongly  coupled  to 
the  lattice. 

Using  a  very  different  approach,  Nguyen,  Nguyen,  and 
Grider3  investigated  the  location  of  traps  in  an  AlGaN/GaN 
heterostructure  FET  by  observing  a  reduction  in  current 
compression  with  increasing  negative  gate  bias.  They  con¬ 
cluded  from  this  that  the  traps  resided  either  in  the  AlGaN 
layer  or  at  the  interface.  A  similar  conclusion  was  deter¬ 
mined  for  AlGaAs/GaAs  modulation-doped  FETs10  over  a 
decade  ago.  Recently,  Trassaert  et  al.4  have  observed  the  ef¬ 
fects  of  surface  states  in  their  GaN  MESFETs.  This  conclu¬ 
sion  was  based  on  the  lack  of  any  difference  in  threshold 
voltage  between  static  and  pulsed  measurement  conditions. 

Because  of  the  simplicity  of  the  MESFET  structure,  the 
traps  in  the  current  investigation  can  reside  only  in  the  HR 
GaN  layer,  at  the  channel/HR-GaN  interface,  or  as  surface 
states.  We  do  not  believe  that  surface  states  or  surface 
charges  can  account  for  the  observed  current  collapse.  First, 
the  gate  lag  response,  which  is  generally  associated  with 
surfaces,27  shows  little  dependence  on  the  magnitude  of  the 
applied  drain-source  voltage  in  GaN  high  electron  mobility 
transistors  (HEMTs).28  This  is  unlike  the  current  collapse 
mechanism  reported  here,  that  requires  a  high  drain-source 
voltage.  Furthermore,  silicon  nitride  passivation  of  GaN 
HEMTs  has  resulted  in  a  change  of  the  gate  lag  response,29 
with  no  effect  on  the  observed  current  collapse  behavior.  In 
addition,  devices  fabricated  on  conductive  buffer  layers  do 
not  exhibit  current  collapse.29  This  is  a  further  indication  that 
this  phenomenon  is  due  to  charge  trapping  in  the  buffer 
layer. 

The  responsible  traps  cannot  reside  in  the  active  region, 
as  empty  electron  traps  would  simply  be  compensated  by  the 
channel  carriers,  with  or  without  the  high  drain-source  volt¬ 
age.  Their  effect  would  be  only  to  diminish  the  carrier  con¬ 
centration  in  the  channel.  In  this  work,  we  have  assumed  that 
the  traps  are  in  the  bulk  of  the  HR  GaN.  Such  traps  can  result 
from  point  defects  within  single  GaN  grains,  or  they  can  be 
associated  with  extended  defects,  such  as  deep  centers  local¬ 
ized  within  the  core  of  a  threading  edge  dislocation.  As  we 
noted  earlier,  both  of  the  traps  that  were  observed  through 
photoionization  spectroscopy  exhibit  the  same  optical  ab¬ 
sorption  signatures  as  those  of  defects  observed  to  be  respon- 
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sible  for  PPC,  and  it  is  very  likely  that  they  are  the  same 
defects.  In  both  cases  the  PPC  studies  were  carried  out  on 
thick  GaN  layers  grown  on  sapphire.1314  No  conducting 
channel/HR  GaN  interface  was  present.  It  is  therefore  un¬ 
likely  that  the  traps  studied  in  this  work,  which  we  believe  to 
be  the  same  PPC  centers,  are  located  at  the  interface.  If  it 
were  assumed,  however,  that  the  traps  reside  only  at  the  in¬ 
terface,  it  is  straightforward  to  show  that  the  analysis  de¬ 
scribed  in  Secs.  IV  and  V  remains  essentially  unchanged: 
The  areal  charge  density  eN~(t)xT  in  the  charged  region  of 
the  HR  GaN  is  replaced  by  eM{t),  where  Af(t)  is  the  areal 
density  of  filled  traps  at  the  interface. 

VII.  SUMMARY 

The  dependence  upon  light  intensity  (exposure)  of  the 
reversal  of  current  collapse  was  measured  in  a  GaN 
MESFET  for  two  distinct  trapping  centers.  In  order  to  extract 
meaningful  parameters  from  the  data,  the  process  was  mod¬ 
eled  as  resulting  from  the  photoionization  of  trapped  carriers 
in  the  HR  GaN  layer  and  the  subsequent  return  of  these 
carriers  to  the  channel  under  the  influence  of  the  built-in 
electric  field  at  the  HR/channel  interface.  The  modeling  pro¬ 
cedure  also  determined  the  conditions  under  which  the  opti¬ 
cal  response  function  S{hv),  that  measures  the  reversal  of 
current  collapse,  is  proportional  to  the  photoionization  cross 
section  of  the  trap:  These  conditions  are:  (1)  that  the  device 
is  operated  in  the  linear  regime,  and  (2)  that  the  measure¬ 
ments  are  carried  out  at  low  optical  excitation.  The  effect  of 
the  recapture  of  the  photoionized  carriers  on  the  validity  of 
the  proposed  model  was  also  considered.  It  was  found  that 
for  the  traps  considered  here,  which  are  expected  to  have  a 
relatively  small  capture  cross  section,  carrier  capture  be¬ 
comes  important  only  at  long  illumination  times,  when  vir¬ 
tually  all  of  the  trapped  carriers  have  been  returned  to  the 
channel.  Using  the  model  developed  here,  the  experimental 
data  were  fitted,  yielding  a  capture  cross  section  and  an  areal 
trap  concentration  for  each  of  the  traps:  crl  =  6X  1 0  17  cm2, 
Af^XlO11  cm'2;  <t2£=2.8X  10'15  cm2,  U2=  6X1011 
cm'2.  Thus,  the  two  traps  occur  with  comparable  concentra¬ 
tions,  and  contribute  comparably  to  the  current  collapse.  The 
large  difference  in  their  measured  optical  response,  S{liv), 
was  found  to  be  due  to  the  large  difference  in  their  photo¬ 
ionization  cross  sections. 
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